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ABSTRACT 

Using the results of a numerical simulation which follows the evolution, metal en- 
richment and energy deposition of both Population III and Population II stars, we 
predict the redshift dependence of the formation rate of black hole remnants of Popu- 
lation III stars with masses 100 — 500Mq and of neutron stars (black holes) remnants 
of Population II stars with masses 8 - 2OM (20 - 4OM ). 

We describe the gravitational wave spectrum produced by Population III and 
Population II sources adopting the most appropriate signals available in the literature 
and we compute the stochastic backgrounds resulting from the cumulative emission 
of these sources throughout the history of the Universe. 

With the aim of assessing whether these backgrounds might act as foregrounds 
for signals generated in the Inflationary epoch, we compare their amplitudes with the 
sensitivity of currently planned and future ground/space-based interferometers. 

The predicted Population III background lies in the sensitivity range of Ultimate- 
DECIGO, adding as a confusion-limited noise, with a peak amplitude of flcwh 2 — 
3 x 10~ 15 at / = 2.74 Hz. However, differently to previous claims, we find that the 
gravitational wave background generated in the Inflationary epoch may dominate for 
/ ^ 2 Hz. At frequencies / ^ 10 Hz, the background generated by Population II 
stellar remnants is much larger than that associated with Population III stars, with 
peak amplitudes ranging between 10 -12 ^ flowh 2 ^ 7 x 10~ 10 at frequencies / G 
(387 — 850) Hz: progenitors with masses in the range 20 — IOOMq contribute with a 
nearly monotonically increasing behavior, whereas for stars with masses 8 — 25Mq the 
resulting background shape depends on the waveforms used to represent the collapse 
signal. 

Finally, we explore a scenario in which Super Massive Stars, formed in proto- 
galactic dark matter halos out of gas of primordial composition and irradiated by a 
strong ultraviolet background, collapse to Super Massive Black Holes (SMBHs). Even 
assuming, as an upper limit to their formation rate, that the mass density of these 
SMBHs at z = 10 equals the presently observed value, the resulting gravitational wave 
background is too low to be detected from the space-based interferometer LISA; it is in 
the sensitivity range of Ultimate-DECIGO but in a frequency region seriously limited 
by the Galactic binary confusion background. 

Key words: gravitational waves - galaxies: formation -stars: early type - cosmology: 
theory. 



1 INTRODUCTION 

Gravitational waves of cosmological origin could be the re- 
sult of a large variety of astrophysical and cosmological pro- 
cesses that develop in the very early Universe. As a con- 
sequence, our high redshift Universe is expected to be per- 
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meated with a background of gravitational radiation. De- 
pending on their origin, these stochastic gravitational wave 
backgrounds (GWBs) will show different spectral properties 
and features that it is important to investigate in view of a 
possible, future detection. 

At present, several gravitational waves interferometers 
(VIRGO, LIGO, GEO600, TAMA) are operating and tak- 
ing data in the frequency range ft! 10 Hz — 3 kHz. In the 
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near future, the low-frequency window will be accessible 
to observations: the Laser Interferometer Space AntenncQ 
(LISA) covering the frequency range (O.l-lOO)mHz will be 
launched, and next-generation space detectors targeting 
(O.l-l)Hz as Decihertz Interferometer Gravitational wave 
Observatory (D ECIGO) and Big Bang Observ ed (BBO) are 
being planned (|Setoll2006l ; [Kudoh et alj|2006h . 

These experiments will have as possible objectives to 
detect or to set bounds on primordial GWBs, originating 
at the Big Bang, as well as backgrounds produced by the 
cumulative emission of astrophysical sources throughout the 
history of the Universe. In particular, the low-frequency re- 
gion is where we expect the contribution of the GWB from 
Population III stars (Pop III) that are the first stars to form 
in the Universe. Pop III stars do not contain metals and 
their formation and evolution may be different from th at of 
stars of late r generations (see iBromm fe Larson! 12004 and 
lGloverll2005l for thorough reviews of the subject). In partic- 
ular, recent theoretical studies suggest that below a criti- 
cal metallicity of Z CI = 10 -5±1 Zq the reduced cooling effi- 
ciency and large accretion rates favor the formatio n of mas- 
sive stars, with characteristic masses > lOOM^ dOmukail 
200d;lBromm et. a|200ll ; ISchneider et~alll2002l . I2003l . l2006al ; 
Qmukai et alj|2005h . 



These stars are predicted to collapse to black holes of 
comparable masses (except in the mass range 140 — 260A/q 
where they are com pletely disrupted in pair -instability su- 
pernova explosions, iHeger fc Wooslevl l2002h and therefore 
are expected to be efficient sources of gravitational waves. 

An earlie r estimate of the GWB produc ed by these stars 
was d one by ISchneider et alj l|2000l) and Ide Arauio et al.l 
<|2002l ). in addition to that a ss ociated to core-collapse super- 
novae (|Ferrari et all Il999al lbl; Ide Arauio et al.ll2004) . How- 
ever, these earlier studies were limited by our poor under- 
standing of the characteristic ma sses and formation rate s 
of Pop III st a rs. M ore r ecently, Buonanno et all (2005), 
ISandick et ail |2006h and ISuwa et al.l (l2007al ) pointed out 
that in some regions of the parameter-space and depend- 
ing on the adopted Pop III star formation rate, the GWB 
produced by the collapse of these first stars could be com- 
parable or mask almost completely the primordial back- 
ground predjcted_bj^^andArd Inflatio nary models. In ad - 
dition, iBuonanno et al.1 (120051 ) and ISandick et ail (|2006l ) 
have taken into account the contribution associated to core- 
collapse supernova events (the end-product of Pop II stars 
with progenitor masses > 8Mq) in order to investigate in 
which frequency ranges these two components could be dis- 
criminated. 

In order to give a realistic estimate of the GWB spec- 
trum produced by the formation of the black holes remnants 
of the first stars, it is necessary to (i) reduce the big uncer- 
tainty related to the prediction of Pop III star formation rate 
and (ii) adopt suitable models for the gravitational wave sig- 
nals emitted during the collapse of these primordial stars. 

This paper is an attempt to satisfy these two re- 
quirements. We compute the GWB produced by rem- 
nants of Pop III stars in the mass range ( 100-500) Mq us- 
ing the Pop III cosmic star formation rate obtained by 
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iTornatore et al.l l|2007al ) in a recent numerical simulation 
which follows the evolution, metal enrichment and energy 
deposition of both Pop III and Pop II stars, and give a reli- 
able estimate of their relative rates at high redshift. As will 
be discussed later, we adopt as a template for the gravita- 
tional waveforms associated to the collapse of P op III stars 
the res ults of recent numerical simulations by ISuwa et al.l 
(2007a) and discuss the relate d uncertainties comparing 
with other independent studies (jFrver et al . 2001; L iu et al.l 
120071 : ISuwa et al.ir2007bl ) . We also evaluate the contribution 
to the GWB of Pop II stars, which explode as core-collapse 
supernovae leaving a neutron star (progenitor masses in 
the range [8 — 20] Mq) or a black hole remnant (progen- 
itor masses between [20 — 40]A/q), or which directly col- 
lapse to black holes (progenitor masses larger than 40Mq). 
In these ranges of Pop II progenitor masses we choose the 
most appropriate gravitational wave signals available in the 
literature. 

The plan of the paper is as follows. In Section 
[2] we bri efly describe the num erical simulation per- 
formed by ITornatore et al.l (|2007aT ) and the model for the 
Pop Ill/Pop II cosmic star formation rate evolution that we 
have adopted in our analysis. In Sections [3] and [4] we sketch 
out the main features of the waveforms which describe the 
gravitational emission of single Pop III and Pop II sources. 
In Section [5] we present the density parameter, Hgw, of the 
GWBs, discuss their detectability by existing and planned 
detectors and the possibility that they might be foregrounds 
limiting the detection of primordial GWBs. In Section[6]we 
explore the GWB produced by Super Massive Star formed 
in metal- free protogalaxies at redshift z = 10. Finally, in 
Section [7] we draw our conclusions. 

Throughout our work we have adopted a ACDM cos- 
mological model with parameters !Jm = 0.26, Qa = 0.74, 
h = 0.73, Qb = 0.041, in agreem ent with the three-year 
WMAP results (|Spergel et al.ll2007l ). 



2 THE POPULATION HI/POPULATION II 
COSMIC STAR FORMATION RATE 

Our present understanding of Pop III star formation sug- 
gests that these stars are not necessarily confined to form 
in the first dark matter halos at redshift z > 20 but may 
continue to form during cosmic evolution in regions of suffi- 
ciently low metallicity, with Z < Z CI . Thus, the Pop III star 
formation rate and the cosmic transition between Pop III 
and Pop II stars is regulated by the rate at which metals 
are formed and mixed in the gas surrounding the first star 
forming regions, a mechanism that we generally refer to as 
chemical feedback. 

Semi-analytic studies which implement chemical feed- 
back generally find that, due to inhomogeneous metal en- 
richment, the transition is extended in time, with coeval 
epochs of Pop III and Pop II star formation, and that 
Pop III stars c an continue to form down to moderate red- 
shifts, z < 5 ( [Scannapicco et al l 120031 : iFurlanetto fc Loebl 
120051 : ISchneider et al.ll2006bl ). Yet, a direct observational ev- 
idence for the existence of Pop III stars is still lacking. 

To better assess the val idity of these semi-analytic mod- 
els, ITornatore et al.l (|2007al ) have performed a set of cosmo- 
logical hydrodynamic simulations with an improved treat- 
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ment of chemical enrichment (|Tornatore et al,ll2007tj ) . In the 
simulation, it is possible to assign a metallicity-dependent 
stellar initial mass function (IMF). When Z > Z cv Pop II 
stars are assumed to form according to a Salpeter IMF 
$(M) tx M~ < - 1+x) with x = 1.35 and lower (upper) 
mass limit of O.IMq (IOOA/q). Stars with masses in the 
8 — 40Mq range explode as Type-II SNe and the explo- 
sion energy and metallici t y-dep endent metal yields are taken 
from lWooslev fc Weaver! jl995l ). Stars with masses > 40A/q 
do not contribute to metal enrichment as they are assumed 
to directly collapse to black holes. Very massive Pop III 
stars form in regions where Z < Z CT ; since theoretical mod- 
els do not y et provide any indicatio n on the shape of the 
Pop III IMF. lTornatore et all |2007a|) adopt a Salpeter IMF 
shifted to the mass range 100 — 500Mq; only stars in the 
pair- instability range (140 — 260Mq) contribute to metal- 
enrichment and the metal yie l ds an d explosion energies are 
taken from iHeeer fc Wooslevl (|2002f) . The simulation allows 
to follow metal enrichment properly accounting for the fi- 
nite stellar lifetimes of stars of different masses, the change 
of the stellar IMF and metal yields. The numerical schemes 
adopted to simulate metal trans port ad diffusion are dis- 
cussed in iTornatore et~aH (|2007bh to which we refer the in- 
terested reader for more details. 

The top panel in Fig. [1] shows the redshift evolution 
of the cosmic star formation rate and the contribution of 
Pop III stars predicted by the simulatior0. In this model, the 
critical metallicity which defines the Pop Ill/Pop II transi- 
tion is taken to be Z CI = 10~ 4 Zq, at the upper limit of the 
allowed range. However, additional runs show that decreas- 
ing the critical metallicity to W~ 6 Zq reduces the Pop III 
star formation rate by a factor < 10. It is clear from the 
figure that Pop II stars always dominate the cosmic star 
formation rate; however, in agreement with previous semi- 
analytic studies, the simulation shows that Pop III stars 
continue to form down to z < 5, although with a decreasing 
rate. Over cosmic history, the fraction of baryons processed 
by Pop III stars is predicted to be ft — 2 x 10 -6 . 

The number of gravitational wave sources formed per 
unit time out to a given redshift z can be computed inte- 
grating the cosmic star formation rate density, p*(z), on the 
comoving volume element and restricting the integral over 
the stellar IMF in the proper range of progenitor masses, 
that is 
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where the factor (1 + z) at the denominator takes into ac- 
count the time-dilation effect and the comoving volume ele- 
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Figure 1. Top panel: redshift evolution of the comoving star 
formation rate density (solid lines). The dotted line shows the 
contribution of Pop III stars (see text). Bottom panel: redshift 
evolution of the number of gravitational wave sources formed per 
unit time within a comoving volume. The top curves show the 
contribution of Pop II stars leaving behind neutron stars (NS,II 
solid line) and black holes (BH,II triangles); the bottom dotted 
line represents the contribution of Pop III stars collapsing to black 
holes (BH,III see text). 



ment can be expressed as, 



dV = Amr I —J e(z)dz 
e(z) = [Q M {1 + zf + Q A ] 



(4) 



Because of the different gravitational wave emission mech- 
anisms, we have explicitly divided Pop II sources between 
those that lead to the formation of a neutron star remnant 
(NS,II) and those which collapse to black holes (BH,II). 
Pop III sources have been considered to be all stars which 
collapse directly to black holes (BH,III). Therefore we cut 
the mass range 140Mq ^ M sj 260Mq of pair instability 
supernovae because in this range of masses the progenitors 
encounter the electron-positron pair creat ion instability and 
are co mpletely disrupted by the explosion (|Heger fc Wooslevl 
2002). Note, however, that these results strictly apply only 
to non-rotating stellar models and that rotation may af- 
fect mass-loss, nucleosynthesis, and the pre-supernova struc- 
ture. Recent analysis show that, at least for metal- free stars 
(Z = 0), up to 10% of the initial stellar mass can be lost 
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when the critical break-up limit is reac hed but this is not im - 
portant enough to change their fate (jEkstrom et al.ll200Sh . 
The results are shown in the bottom panel of Fig. [1] It is ev- 
ident that the rate of Pop III stars collapse to black holes at 
z 6 is ~ 10 -5 times smaller than the rate of gravitational 
wave sources associated to Pop II stars. The latter are dom- 
inated by core-collapse SNe leaving neutron star remnants, 
because of the adopted shape of the stellar IMF. 

To compute the gravitational signal emitted by each 
single source, it is important to estimate the fraction of the 
initial stellar mass which ends up in the final remnant NS or 
BH. For Pop III stars, we will neglect mass loss and assume 
that the final BHs have the same mass of their progenitor 
stars (100 — 500Mq). For Pop II stars with masses in the 
Type-II SNe progenitor range (8 — 40Mq), the mass of the 
remnant is computed according to the original Woosley & 
Weaver (1995) grid; for larger masses (> 40Mq) the mass 
of the final BH is taken to be equal to the helium core mass 
before the collapse, M Hc = ||(M - 20M Q ). 
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Figure 2. Spectral distribution of GWs calculated inlSu wa et all 
j2007al) (model S07, solid line) and in iBuonanno et al.l (120051) 
(model B05 with E v = 10 55 erg and < q >= 0.03, dotted line). 



3 POP III STARS AS GRAVITATIONAL WAVE 
SOURCES 

Recent studies have shown that the collapse of Pop III 
stars to black holes could be a much more efficient source 
of gravitational w a ves than today's su pernovae populations 
l|Frver et all l200ll ; ISuwa et all l2007al ). An estimate of the 
total energy released in gravitational w aves during the firs t 
few seconds of collapse has been done in IFrver et al.1 l|200lf ). 
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In their two-dimensional numerical simulations the authors 
have considered the entire life and death of a zero-metallicity 
star of 300Mq, taking the effects of general relativity 
and neutrino tran sport into account. As first noted in 
IFrver et al l (|200lD , the main features of the Pop III col- 
lapse are different from those that appear in ordinary core- 
collapse SNe. In the latter case, core collapse is halted when 
the central density exceeds ~ 10 14 g cm~ 3 and when neutron 
degeneracy and nuclear forces become relevant. Conversely, 
the core of a 300Mq star is so large that it collapses into 
a black hole before nucl ear forces can affect it (jSuwa et al.l 
l2007al ; IFrver et al.ll200ll ). 

In add i tion, as noted by IFrver et al.l l|200ll ) and 
ISuwa et al l (l2007al l. rotation and higher central tempera- 
tures in Pop III progenitors, may halt the collapse of the 
core and produce a weak thermal bounce at lower central 
density of a few times ~ 10 12 g cm -3 . After this bounce, 
the final fate of hot neutronized core, as soon as entropy 
gradients becomes negligible, is to collapse to a black hole. 

Lacking mo re detailed numer i cal mo dels on the collapse 
of Pop III stars, IBuonanno et alj (2005) have modelled the 
gravitational wave signal emitted in the collapse of a 300Mq 
Pop III star using as a template the waveform of an ordinary 
core-collapse supernova, 
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as obtained bv lMuiier et~ai1 (|2004h using sophisticated two- 
dimensional numerical simulations. The stellar core-collapse 
is taken to be axially symmetric and is based on a detailed 
imple mentation of Boltz mann's equation for neutrino trans- 
port (Mul ler et all [2004). In equation (©, h is the Fourier 
transform of h, G is Newton's constant, E v is the total en- 
ergy emitted in neutrinos, < q > is the average value of the 
anisotropy parameter (defined in lMuller fc Ja nka 199?]), r is 
the distance of the supernova event, a and b are respectively 
200Hz and 300Hz. This choice of a and b values is motivated 
in order to r eproduce the spectr al shape of the sl5r model 
simulated in lMuller et~aH (|2004T l. 

The energy spectrum adopted by IBuonanno et alj 

l|2005h therefore is 



dE G w 16i7 2 c 3 r 2 2,wn,2 



[erg/Hz] 



(7) 



with the additional requirement that the total ener gy emit- 
ted had the same value found bv lFrver et al.l (|200ll ). 



Egw = 



15G 



dff\h(f)\ 2 = 2xlO- 3 M c 2 . (8) 



The resulting GW spectrum is shown in Fig. [2] 

More recently, ISuwa et al.1 (|2007al ) have studied the ro- 
tational collapse of Pop III stars in the same mass range. 
These studies have been performed using two-dimensional 
hydrodynamic simulations, using a realistic equation of 
state based on the relativistic mean field theory ( Shcn ct al. 
1998), and taking neutrino cooling and state-of-the-art re- 
actions of neutrin os into account. F or comparison with the 
previous study of IFrver et al.l (|200ll ). they discuss the GW 
emission from the collapse of a 300AfQ Pop III star and 
assume the same initial core configuration. The total en- 
ergy emitted in gravitational waves is found to be, Egw ~ 
0. 2 x 10~ 3 M^c 2 , th at is a factor 10 smaller than the result 
of IFrver "et al. (2001) (see eq(5}, and the authors ascribe this 
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discrepancy to the assumed different initial angular momen- 
tum distributions. 

In Fig. [2] the solid line shows the GW spectral distri- 
bution obtained in ISuwa et al l (|2007ah . As it is clear from 
the comparison of the two curves, the contribution to GW 
emission from neutrinos in the model of ISuwa et al.l (|2007al ) 
dominates the spectrum below 10 Hz and the matter contri- 
bution peaks at around 20 Hz, whereas adopting the spec- 
tral energy distr i bution of ordinary supernovae, as done by 
iBuonanno et al.l (|2005T l. these frequencies are shifted in the 
kHz range. 

As discussed by iFrver et al.l (|200lf ) and ISuwa et al.l 
l|2007al ). these differences are expected, given the different 
nature of the progenitor and the distin ct collapse dynamics . 
It is interesting to note, as shown by ISuwa et al.l (|2007bh . 
that since most of the inner core is absorbed into the black 
hole, only about 10% of the gravitational energy of the core 
can be carried away by neutrinos (even in the case of the 
most rapidly rotating model). On the other hand, neutrinos 
carry away 99% of the gravitational energy of proto-neutron 
stars in ordinary core-collapse supernovae. 

Unless otherwise specified, in our estimates of the grav- 
itational wave background produced by Pop III stars, we 
will describe the single source energy spectrum using equa- 
tion Q with f\h(f)\ taken from the results of ISuwa et al.l 
(|2007ah . 



4 POP II STARS AS GRAVITATIONAL WAVE 
SOURCES 

In spite of the many elaborate studies carried out so far, the 
physics involved in core-collapse supernovae is not yet com- 
pletely understood. During the past years, theorists have 
done many efforts to give reliable predictions of the gravita- 
tional wave signals associated to these catastroph ic ev ents 
(for r ecent reviews we refer to iKotake et al. I l2006l andlottl 
l200Sh . 

In order to estimate the stochastic gravitational wave 
background from Pop II stars, we have adopted different 
single source gravitational wave spectra obtained from the 
most sophisticated numerical simulations available in the 
literature. 

The most recent calculations incorporate general rel- 
ativity as well as the most relevant physics, such as a 
microphysical finite-temperature nuclear equation of state, 
a scheme to account for electron capture and neutrino 
losses during the collapse, an accurat e treatment of t he 
neutrino transport ([Burrows et al"1l2006l ; lOtt et aT I l2007al lbl; 
iDimmelmeier et al.|[2007i . l2008l '). 

These studies seem to indicate that, in spite of the 
dependence of the gravitational wave burst on the prec- 
ollapse central angular velocity, on the progenitor mass, 
on the equation of state, etc, all the simulated waveforms 
exhibit a generic bo unce dynamics and signal morp hology 
i|Ott et al.l r2007al lbt IDimmelmeier et al.l |2007| . I2008D . The 
bounce dynamics seems to be governed by the stiffening of 
the equation of state at nuclear density and, as a conse- 
quence, the resulting waveforms show one pronounced large 
spike at the bounce and a gradually damped ring-down 
l|Dimmelmeier et aLlkoOST ). 



4.1 Pop II progenitors with masses in the range 
[8-20] Mq 

To model the single source gravitational spectrum for Pop II 
stars with masses in the range [8-20] Mq , we have used equa- 
tion (|7J) a dopting the f \ h(f) | predicted by the model labelled 
as sl5r in iMuller et al1 j2004). This choice provides the most 
optimistic co nditions for the tot al energy released in gravita- 
tional waves (|Muller et al.ll2004h . and has been considered to 
be representati ve of the emission in co r e-collapse SNe by pre - 
vious studies l|Buonanno et al] 120051 ; ISandick et al.l l2006). 
We have assumed that the total energy emitted in gravita- 
tional waves is Egw — 1-8 x 10 _8 Mqc 2 (that corresponds 
to a conversion efficienc y of the mass energ y in gravitational 
waves of e ~ 1.5 x 10" 



1 01 the mass energy m gravitational 
Muller et al]|2004). The lat est stud- 



ies carried out by IDimmelmeier et all (2007. 200jl) confirm 
that this choice for the single source spectrum is represen- 
tative of core-collapse supernovae with progenitors in the 
above mentioned mass range. 

Recently. iBurrows et al.l |2006) have proposed an acous- 
tic mechanism for core-collapse supernovae associated to 
stellar progenitors in the mass range 8 — 25A/q. In this 
model, g modes (mainly 1=1,2 core g mode) are excited since 
the proto-neutron star undergoes the standing-accretion- 
shock instability (SASI, see for a review on the subject 
IBurrows et "all l2006l : lOtt et al.l [20061 : IOttll200Sft induced by 
turbolence and accretion downstream of an unstable and de- 
formed stalled supernova shock. This mechanism has been 
suggested to be sufficient to drive a supernova explosion. 

InlOtt" et al.l (|2006l ) the authors investigate this mecha- 
nism using a two dimensional axisymmetric newtonian sim- 
ulation. The results seem to indicate that the emission pro- 
cess in core-collapse supernovae may be dominated by the 
oscillations of the proto-neutron star core. The authors ex- 
plore different type of progenitors and find that the more 
massive is the progenitors iron core, the higher are the fre- 
quency and oscillation amplitudes and, as a consequence, 
the stronger is the gravitational wave emission. It is in- 
teresting to note that for these models the contribution to 
gravitational waves of anisotropic neutrino emission is com- 
pletely negligible compared to that associated to core oscil- 
lations. In particular their estimate of the contrib ution of 
the anisotropic neutrino em ission (calculated as in lEpsteinl 
ll978l ; lMuller fc Jankalll997f) is fiv e time s smaller than previ- 
ous one obtained in lMuller et al.l (|2004r ). The authors opin- 
ion is that this is due to the different approach in treating 
neutrino transport. We refer the interested reader to the 
original paper for further details on the subject. 

The total energies emitted in gravitational waves are 
around Egw ~ 10~ 8 Mqc 2 (for a non rotating progenitor 
of 11 A^q and a slowly rotating progenitor of 15 Mq) with 
the exception of a non rotating progenitor of 25 Mq, with a 
very massive iron core (1.92Mq), where the emitted energy 
is as large as Egw ~ 8.2 x 10 _5 Mqc 2 . It is important to 
note that these results are based on 2D Newtonian gravity 
simulations; rotation and 3D simulations may change the 
structure of the SASI and consequently t he energy e mitte d 
in gravitational waves, as point out by lOtt et all {2006). 
As an extreme and promising possibility, we also explore 
the gravitational background produced by core-collapse SNe 
through this mechanism. 

The energy spectra predicted by Muller et alj (|2004h 
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Figure 3. The GW ene rgy power spe c tra, fo r neutron star core- 
collapse, obtained from iMuller et ail 112004 ) ( model M04, with 
E v = 3 X 10 53 erg and < q >= 0.0045) and lOtt et ail d2006h 
(model O06). 
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Figure 4. The GW energy power spectra for m odels D5a, D5d 
and A5b obtained in lSekiguchi k Shibatal j2005l ) (models SS05). 



and lOtt et ail i|2006l ) are shown in Fig. If we compare 
th e two spec t ra, w e see that the energy spectrum obtained 
in lOtt et ail (2006) shows a pronounced gravitational wave 
emission in the frequency range 600 — 1000Hz, the first burst 
of gravitational emission is centered about 800Hz, and the 
second one near 950Hz, due to the excitement of an 1=2 core 
g-mode. 

4.2 Pop II progenitors with masses in the range 
[20-100] M & 

Progenitors with masses between [20-100]Mq may lead to a 
prompt collapse to black hole or to the formation of a proto- 
neutron star that, due to the subsequent fallbac k will end 
its life as a bl ack hole (ISekiguchi fc Shibata|[2005l ). 

Recently, ISekiguchi fc Shibatal (^OOpF have run numer- 
ical simulations of the collapse of rotating massive cores 
associated to stellar progenitor masses in the range [51- 
98]A#q. These simulations have been performed in full 
general relativity using a parametric equation of state 
l|Dimmelmeier et ail 12002a b) that mimic a realistic one 



0. Gravitational waveforms are computed in terms of the 
quadrupole formula taking into account only matter motions 
and neglecting neutrino contribution (we refer the interested 
reader to the original paper for further details). 

For the purpose of the present analysis, we adopt the 
gr avitational wave spect r a from models D5a, D5d and A5b 
of ISekiguchi fc Shibatal i|2005l ). which show the strongest 
gravitational wave emission with efficiencies ranging be- 
tween 0.2 x 10 -6 < e ^ 0.3 x 10" 6 . Despite the fact that 
the selected models are, in some cases, associated to a spe- 
cific choice of the progenitor stellar mass, we will neglect 
this dependence and assume that the same model holds for 
all progenitors in the 20 — IOOMq range. This approxima- 
tion may be less accurate for black holes forming through 
fallback from progenitors with masses < 40Mq; in fact, the 
evolution of the carbon/oxygen envelope as well as neutrino 
cooling (which will play a role for longterm fallback with 
duration > 100 ms) are n ot considered in the simulations 
i|Sekiguchi fc Shibatall2005l ). 

In Fig. |3] we have plotted the energy power spectra for 
models D5a, D5d, A5b. For all three models, the spectra 
peak at / ~ 1 — 2kHz. The energy power spectrum of model 
D5a is the result of a strong bipolar explosion (this feat ure 
is found for the first time in ISekiguchi fc Shibatal [20051 ); it 
is broad in the low frequency region and flatter than those 
of the other two models. All spectra peak at / ~ 1 — 2kHz. 
Due to the smaller progenitor mass (51Mq) the energy spec- 
trum of model A5b shows a lower peak then in model D5d 
(progenitor mass 75Mq) and shifted to lower frequencies. 



5 GRAVITATIONAL WAVE BACKGROUND 
FROM POP III/POP II STARS 

In this section we discuss the contribution of the different 
gravit ational wave sources to the background radiation. Fol- 
lowing [Schneiderletjiil (|2000f ) . the spectral energy density of 
the stochastic background can be written as 



dE 



dE 



dSdfdt J J M . v ' ' x dSdf 
where dR(M, z) is the differential source formation rate 
p*(z) dV 



0) 



dR(M, z) 



(1 + z) dz 



$(M)dMdz, 



(10) 



and ( ) , is the locally measured average energy flux emit- 
ted by a source at distance r. For sources at redshift z it 
becomes, 



( JE_ )= {l + zf dE hw 



(ii) 



where / = / e (l + z)^ 1 is the redshifted emission fre- 
quency / e , and di J (z) is the luminosity distance to the 
source. It is customary to describe the GWB by a dimen- 
sionless quantity, the closure energy density Ogw(/) = 



5 Th is realistic equation of state is the EOS Shen l lShen et al 
Il998t) and we note that this is the sa me EOS used in~ lott*^tTl 
J2006h and lDimmelmeier et ail ll200d) . 
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Figure 5. The function Cgw evaluated by assuming a frequency 
cut off / ^1 Hz. 
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Figure 6. The function SJqw in the zero-frequency limit. 



p cr 1 (dp gm /dlogf), which is related to the spectral energy 
density by the equation 



«GW(/) = 



dE 



dSdfdt 



(12) 



where p cr = 3 H / 8nG is the cosmic critical density. 

In Fig. [5] we show the function SIgw evaluated by as- 
suming a frequency cut-off / ^1 Hz in the single source 
spectrum. This choice is motivated by the fact that current 
numerical simulations describe the source collapse at most 
for a few seconds, and cannot predict the emission below a 
fraction of Hz. We see that below ~ 10Hz the Pop III GWB 
dominates over the Pop II background (except for a very 
small region around ~ 1Hz, where the Pop III GWB has a 
minimum); its maximum amplitude is JIgw^ 2 — 3 x 10~ 15 
at / = 2.74 Hz. At larger frequencies the background pro- 
duced by Pop II stars is much larger than that of Pop III. 
Stars with progenitors in the range (20-100)Mq contribute 
with a nearly monotonically increasing behavior, reaching 
amplitudes flcwh 2 ~ 4 x 10" 10 at / SE (759 - 850) Hz. For 
stars with progenitors in the range (8-25) Mq the GWB has 
a shape which depends on the waveform used as represen- 
tative of the collap se: that obtained using the waveform of 
iMuller et alj lj2004h is quite smooth and peaks at / = 387 
Hz, with amplitude QGwh 2 ~ 10 -12 . The background ob- 
tai ned using wavefo rms produced by the core-collapse model 
of lOtt et al.l f2006), instead, is comparable in amplitude 
with those produced by more massive progenitors, reach- 
ing the maximum amplitude S7gw^ 2 — 7 x 1CP 10 at / = 485 
Hz. 

To evaluate the background in the low frequency re- 
gion / ^ 0.1 Hz, w e can app l y the so called zero-frequency 
limit dSmarrl 1 19771; lEpsteinl ll978l ; IMuller fc Jankal 1 19971 ; 
iBuonanno et al.ll2005h ~ and extend the single source wave- 
form f\h(f) \ to lower frequencies (using equation|S]for model 
M04, and tabulated values for models S07 and O06), where 
the emission is dominated by the neutrino signal (see Fig.0. 

The results are shown in Fig. [6] where we plot QGV/h 2 
for all models shown in Fig. \5 j with the exception of the 
Pop II signals derived from the 



Sekiguchi fc Shibatal (|2005l ) 



single source spectra because in these computations neutri- 



nos contribution is neglected and the extension would not 
be appropriate. 

In the low frequency region, Qgw x f because f\h(f)\ 
tends to a constant as shown in Fig. [2] In addition, the 
Pop III-SIgw dominates over the Pop II, because the GW 
emission increases with the stellar progenitor mass (the more 
massive is the progenitor, the larger is the neutrino luminos- 

ity). 

We c an compare o ur res ul ts with those r ecentl y ob- 
tained by ISuwa et al.1 (|2007al ). ISandick et"aH (|2006l ) and 
IBuonanno et al.l (|2005l ) , who adopt different Pop III star for- 
mation rates and/or different single source GW waveforms. 
We find that our predicted Pop HI GWB is com parable to 
the lower limit reported by ISuwa et al.l (|2007al ); they use 
the same single source waveform, but a diffe rent star forma- 
tion rate derived from ISandick et al. I i|2006l) and artificially 
shifted downward by a factor 7000. For this reason, they 
predict a baryon fraction in Pop III stars of 10 -5 , a value 
which is much closer to our numerical result (see Section [2]). 

Our Pop HI background is smaller than t hat f ound 
by ISandick et alJ (I2 OO6) and IBuonanno et al.l (|2005l ). In 
IBuonanno et al. 1 20051 ). the authors make a crude estimate 
of the Pop III star formation rate, assuming that all Pop III 
stars form in a burst at z = 15 with the same mass (300 
Mq), and varying the fraction of baryons which go into Pop 
III stars in the range 5 x 10 -7 — 10 -3 . Comparing with their 
lower limit background, our result shows a peak at lower 
frequency and with smaller amplitude, and a larger low- 
frequency limit. These differences are naturally explained 
by th e adopted sin g le sou rce waveform (see Fig. UJ. Fi- 
nally, ISandick et al.l l|2006l ) show a cumulative GWB from 
Pop III and Pop II stars for a cosmic star formation model 
where the baryon fraction in Pop III stars is 7 x 10~ 2 , much 
higher than tha t predic ted by the numerical simulation of 
iTornatore et al.l l|2007al ) (see Section [2| . It is not surpris- 
ing, therefore, that their signal associated to Pop III stars 
is significantly higher than in our analysis; furthermore, the 
peak amplitude is shifted to frequencies 50 — 100 Hz due 
to the difference in the assumed single source waveform (see 
Fig. (71 . The differences between the Pop II GWBs are less 
pronounced, as the predicted Pop II star formation rates 
are comparable, and the adopted waveform for progenitors 
with masses 8 — 20Mq is the same; still, due to the range 
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Figure 7. The function Hqw is plotted with the sensitivity 
curves of space (LISA TNG, BBO Grand) and ground-based fu- 
ture detectors (LIGOIII). 
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Figure 8. The functio n Ogw is plotted w ith the sensitivity curve 
of Ultimate-DECIGO llKudoh et al.l2006fi . The horizontal dashed 
line is the G WB spectrum produced during slow-roll inflation 
l|Turnerlll997h . 



of efficiency assumed in the present analysis for the col lapse 
of 20 - lOOAfQ progenitors (|Sekiguchi fe Shibatall2005l ). our 
Pop II GW B is ~ 2 times smalle r in amplitude than that 
predicted bv lSandick et all l|2006l ). 



5.1 Detectability 

In this section we discuss the possibility that Pop II and 
Pop III backgrounds might act as a foreground, limiting the 
sensitivity of the detectors to t he primordial GWB produced 
durin g the inflationary epoch (|Kudoh et al.ll2006l ; ISetoll2006l , 
120081 ). 

In Fig. [7J and Fig. [5] we plot Ogw as a function of 
frequency and the sensitivity curves of future d etectors. In 
Fig. [7j we plot the sensitivity of LIGOIII used in lBuonannol 
(2003) and the detection limit for two proposed missions, 
LISA TNG0 (The-Next-Generation ) and BBO in the so- 
called BBO Grand configuration (|Setd 120061 1 . In Fig. [8] 
we plot the detection limit of DECIGO in the so-called 
Ultimate-DECIGO configuration, whos e sensitivity is onl y 
limited by the standard quantum limit l|Kudoh et al.l l2006). 

All the sensitivity curves are obtained assuming corre- 
lated analysis of the outputs of indepe ndent spacecraft 
(or of different ground-based detectors) (|Buonannd 120031 ; 
iKudoh et ai]|2006l ; ISetdl2006l ), with the exception of LISA 
TNG. It is known that a single detector can only place upper 
limits on the amplitude of a GWB. 

The upper limit of the GWB generated in the Inflation- 
ary epoch, is given by the dashed horizontal line in Fig. [7] 
and [8] It is derived from a recent update of cosmological 
constraints on inflationary models (Cosmic Microwave Back- 
ground measurements taken in combin ation with data from 
the A rcminute Cosmology Bolometer iKinnev et al. ( 200E 



) . The data analysis carried out bv lKinnev et al.l 



2006 



2008) is consistent with a primordial spectrum with no run- 



ning spectral index and zero tensor amplitude. The authors 
place an upper limit to the tensor/scalar ratio r < 0.35 at 
95% confidence level. According to these results, we have 



http:/ /www. srl.caltech.edu/ shane/sensitivity/ 



evaluated S2gw from equation (6) of iTurnerl |l997l ). with 
r = 0.3 and no running of the spectral index of tensor per- 
turbation. 

As shown in Fig. LIGOIII, BBO and LISA TNG 
have no chance to detect either the background produced 
by Pop Ill/Pop II stars and by inflation. From Fig.[S]we see 
that the astrophysical backgrounds are within the detection 
range of Ultimate-DECIGO. However, they are smaller than 
the background generated in the Inflationary epoch, which 
clearly dominates. The Pop Ill/Pop II background becomes 
dominant only for / ^ 3 Hz. It is worth reminding that the 
estimate of the Inflationary GWB shown in Fig. [7J and [5] is 
an upper limit. 



5.1.1 Duty Cycle 

An important parameter which quantifies the efficiency of 
the emitters to generate a continuos signal in time domain 
is the duty cycle, D. This is defined as the ratio between the 
typical duration of the signal emitted by a single source, and 
the average time interval between two successive episodes of 
emission. When D — fl, the overall signal is continuous; con- 
versely, if D <C 1, the resulting background is not continuous 
but rather characterized by a shot-noise structure. 

An interesting feature of the GWBs produced by the 
collapse of Pop Ill/Pop II progenitors is that these back- 
grounds generate a s ignal with a p e culiar s hot-noise charac- 
ter, a s suggested by iFerrari et alj (|l999al ) ; ISchneider et al.l 
(2000). In this case, the detector will receive a stationary se- 
quence of bursts, with typical separations much longer than 
the average duration of each single burst. The duty cycle 
measured by a local observer and generated by all sources 
within the comoving volume out to a redshift z can be writ- 
ten as 



D{z) 



dR(z)AT gw {l + z) 



(13) 



where Ar gw is the average time duration of individual sig- 
nals and dR(z) is the number of sources formed per unit 
time at redshift z. We have calculated the duty cycle for 
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Table 1. Duty Cycle for Pop Ill/Pop II progenitors. 



D 

NS(Pop II) BH(Pop II) BH(Pop III) 
8.56X10" 2 3.10X10" 2 1.57X10" 4 



Pop II and Pop III GWBs, assuming that the typical Ar gw 
for these events is the dynamical timescale at core bounce, 
~ 1 ms and ~ 50 ms respectively. The results are illustrated 
in Tableland show that D ~ 1CT 2 for Pop II GWBs and 
~ 10" 4 for Pop III GWBs. 

This shot noise structure of Pop Ill/Pop II GWBs might 
enhance the ability to discriminate these backgrounds from 
the dominant instrumental noise and from signals generated 
in the early Universe, even in frequency ranges where the 
latter have a higher amplitude. 
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Figure 9. The spectral strain amplitude of SMSs signal (dashed 
line) compared to the sensitivity curve of LISA and Ultimatc- 
D ECIGO. The Galactic binary co nfusion background as reported 
in lSath vaprakash fc Schutj d2009t) is shown by the solid line. 



6 GRAVITATIONAL WAVE BACKGROUND 
FROM SUPER MASSIVE STARS 

Recently, it has been suggested that Super Massive Black 
Holes (SMBH) with masses of » 10 5 - 10 6 M Q may form 
in protogalactic dark matter halos with virial tempera- 
tures T v i r 10 4 K provided that (i) the gas has primor- 
dial composition an d (ii) it is irradiated by a strong ul- 
traviolet (UV) flux |Bromm fc Loebl 120031 ). In these con- 
ditions, 3D hydrodynamics simulations show that, due to 
suppressed H2 formation, the gas cools to temperatures that 
are only somewhat smaller than the virial temperature of 
the host halo and condenses isothermally into very large 
clumps, with no sign of fragmentation. As a result, the pri- 
mordial cloud collapses into a central compact object con- 
taining > 10% of the total baryonic mass, possibly through 
an intermediate phase of Super Massi ve Star (SMS) forma- 
tion, (we refer the interested reader to lBromm fc Loebll2003l 
for further details). This scenario would provide potentially 
detectable sources of gravitational waves and an alterna- 
tive way to explain the presence of SMBHs with masses 
« 10 9 M^ powering bright q uasars at redshifts as high as 
z > 6 (|Volonteri fc R ccs 2006). However, high redshift dark 
matter halos with T v \ r ~ 10 4 K are likely already enriched 
with at least trace amounts of metals and dust produced by 
prior star formation in their progenitors. Indeed, it has been 
shown that, even in the presence of a sufficiently strong UV 
background, gas fragmentation is inevitable above a critical 
metallicity, whose value is between Z CT ~ 3 x 10 _4 Zq (in 
the absence of dust) and as low 10~ 6 Zq for a 

dust-to-gas mass ratio of about 0.01Z/Zq (Omukai, Schnei- 
der & Haiman 2008). Therefore, the actual rate of formation 
of (SMSs) SM BHs remains uncertain . 

Following iBromm fc Loebl (|2003h , we can place an up- 
per limit on the SMBH formation rate at z — 10 requiring 
that their total mas s density does not exceed the present- 
day o bserved value dMarconi et al 1 12004 iMerloni fc Heind 
2008 and references therein), that is 

/ smbh f°° dn{ ^!' W) M smbh (M)dM (14) 

iAf min (z=10) dM 

sS 4.3 x lO 5 (/i/O.7) 2 M Mpc" 3 , 



where / sm bh is the fraction of halos that contain a SMBH, 
A^smbh is the SMBH mass which, following the results of 
the simulations, is assumed to scale with the halo mass as 
A/smbh = Q.1(£Ib/£Im)M % n(M, 10) is the Press-Schechter 
halo mass function at redshift z = 10, and M m i n (z = 10) is 
the dark matter halo mass corresponding to a virial temper- 
ature of T vir = 10 4 K at 2 = 10, 

M min (, = lO) = l O W (4)- 3/2 (^) 1/2 

(i is the molecular weight, A c = 18-7T 2 + 82d — 39d 2 , d = 
Qjf — 1 and 

_ Q M (l + z) 3 

M n M {i + z) s + nA 

for a fiat cosmological model. For the adopted cosmological 
parameters, we find that / sm bh = 0.14, that is about 14% 
of protogalactic halos must verify the conditions suitable 
for SMBH formation to reproduce the presently observed 
mass density. From eq. (|15p we can also estimate the average 
SMBH mass, which is found to be ~ 7.5 x 1O 6 M . The grav- 
itational co l lapse o f a SMS to a SMBH has b een s tudied by 
ISaiio et all i|2002t ). IShibata fc Shapirol (|2002| ) and lLiu et at] 
|2007l ). In these papers, the authors suggest that the col- 
lapse of a SMS to a SMBH may lead to the emission of an 
initial grav itational wave bu rst signal. In lShibata fc Shapirol 
l|2002l ) and lLiu et all (|2007l ) the authors investigate the ax- 
isymmetric collapse of a uniformly rotating SMS near the 
mass-shedding limit at the onset of radial instability. These 
simulations are performed in full general relativity using a 
polytropic equation of state (with T — 4/3) and neglect- 
ing neutrino emission and transport (probably not relevant 
for very massive progenitors bec ause of their lo w density 
and temperature, as suggested bv lLiu e t al. 2007). The au- 
thors estimate that ~ 90% of the total mass of the SMS 
is swal lowed by the SMBH at the end of the collapse. Fol- 
lowing ISchneider etafl (|2000l ). we have assumed that the 
single source spectrum emitted during the collapse of SMS 
can be modelled as a Lorenzian (for more details see also 
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ISchneider e7alll2000f) . 



dE_ 



f N (/-/ ) 2 + r2 



(15) 



where f = c/10R g with R g = 2GM smbh /c 2 , T = (27rAt) _1 , 
At = l//o is the duration of the burst, foN is the normal- 
ization and 



AE G w = eM smbh c 



(16) 



is the gravitational energy released in a burst, 

M mhh = 0.1(n B /n M )M mi n(z = 10) (17) 

and e is the efficiency. We assume that the SMS collapses 
to SMBH of comparable mass with an efficie ncy of e = 
2.0 x 10~ 5 . The adopted e value was found by iFrver et al.l 
(2001) in the collapse scenario of very massive progenitors 
(see Section [3}. This is a rough approximation of the grav- 
itational wave spectrum but, in this uncertain evolutionary 
scenario, gives us the possibility to catch its main features. 
Using the estimated SMBH rate and the single source en- 
ergy spectrum described above, we can compute the gener- 
ated GWB using eq. In Fig. we plot the spectral strain 
amplitude and the sensitivity curves of LISA and Ultimate- 
DECIGO. The estimated signal is too low to be detected by 
LISA and it lies above the foreseen sensitivity of Ultimate- 
DECIGO in the frequency range 0.1 mHz^ / ^1 mHz; how- 
ever, in this frequency interval the sensitivity of Ultimate- 
DECIGO (as well as that of LISA) is limited by the unre- 
solved background produced by the gravitational wave emis- 
sion o f Galactic comp act binaries that acts as a confusion 
noise (lNelemansll2009D. shown in Fig. [5] with a solid line 
l|Sathvaprakash fc Schutzll2009h . A possible future detection 
of this signal with Ultimate-DECIGO would require the ap- 
plication of sophisticated algorithms to data analysis, similar 
to those that have been p roposed for the LISA experiment 
l|Crowder fc Cornish|[2007l ). 



7 DISCUSSION AND CONCLUSIONS 

In this paper, we estimate the GWB produced by the col- 
lapse of Pop Ill/Pop II progenitors. We use a new cosmic 
star formation histor y obtained from a recent numerical sim- 
ulation performed bv lTornatore et all (|2007aT ) . We find that 
our Pop Ill/Pop II GWB is below the sensitivity range of 
space detectors like BBO and LISA TNG in the present 
(proposed) configuration, but is in the sensitivity range 
of Ultimate-DECIGO, adding as a confus ion-limited com- 
ponen t . Differently to prev ious results of iBuonanno et ail 
l|2005l ). ISandicket al.1 (|2006l 'l and lSuwa et al l (|2007aT ) we find 
that for / ^ 2 Hz, Pop III GWB is masked by the GWB 
generated in the Inflationary epoch. 

Clearly, the predicted amplitude of Pop Ill/Pop II 
GWB depends on the adopted star formation model and 
on the assumed single source GW spectrum, and we want 
briefly to point out the related uncertainties. We have 
assumed, as a template for GW emission associated to 
Pop HI stel l ar colla pse, the waveform recently obtained by 
ISuwa et al.l (|2007aT ) . The waveforms they find exhibit fea- 
tures that are significantly different if compared to those 
found in the ordinary core-collapse SNe. As noted by the 
authors, their study is based on a Newtonian simulation, 



whereas a detailed understanding of Pop III stars col- 
lapse would requi re a fully general relativistic approach 
i|Suwaet al.ll2007ah . 

In addition, the assumed Pop II GW spectra are affected 
by the incomplete understanding of the SN explosion mech- 
anism, by the complexity of the physics involved (equation 
of state, structure of the stellar core, different approaches 
in treating neutrino transport, etc.) and by different ap- 
proaches in numerical modeling. 

The variations among different single source gravita- 
tional wave spectra lead to comparable uncertainties in the 
corresponding stochastic backgrounds: as an example, a two 
orders of magnitude difference is found for the peak ampli- 
tude of the signal associated to Pop II stars leading to neu- 
tron star remnants , depending on whether the sp ectrum by 
iMuller et all l|2004 ) or the upper limit provided bv lOtt et al.l 
(j200d ) is adopted. 



It is important to note that while the estimated back- 
grounds depend on the Pop III star formation model, such 
as the adopted IMF and the value of Z CI , the advantage of 
our approach is that it allows a self-consistent description of 
the star formation and chemical evolution. Moreover, varia- 
tions of the adopted Z CI parameter within the allowed range 
(10 -6 — 10 _4 Zq) lead to a factor < 10 decrease of the result- 
ing Pop III star formation rate. As long as Pop III stars are 
assumed to be very massive, with M > IOOMq, the result- 
ing star formation rate is not very sensitive to the adopted 
IMF; this is because only stars with masses in the relatively 
small pair-instability progenitor range (140 - 260 Mq) con- 
tribute to metal enrichment. If Pop III stars are assumed to 
form with masses < IOOMq, the smaller metal yields lead 
to a larger number of Z < Z CI regions and therefore to a 
larger Pop III star formation rate. Finally, additional sim- 
ulations with different box sizes and resolution show that 
the scatter in the Pop H i/Pop II ratio remains always < 1% 
ijTornatore et al.ll2007ah . 

In addition, we have given an upper limit to the GWB 
produced by the collapse of SMSs to SMBHs, events that 
may occur in metal-free halos with virial temperatures 
Tyjr ^ 10 4 K at z = 10 exposed to a strong UV background 
ijBromm fc Loebll2003l ). LISA has no chance to detect this 
low signal and Ultimate-DECIGO, in this frequency region, 
is seriously limited by Galactic binary foregrounds. 

In conclusion, we want to point out that if future space 
missions, like Ultimate-DECIGO, will detect the Inflation- 
ary GWB, it will be necessary to disentangle this back- 
ground from those produced by astrophysical sources. 



Recent stud ies performed bv lDrasco fc Flanagan! (|2003l ) 
and ISetol I 2008) are moving in this direction suggesting dif- 
ferent detect i on me thods for non-Gaussian GWB. In par- 
ticular ISetol (2008) has proposed a method to check if 
the detected GWB is an inflation-type background or if it 
is contaminated by undetectable weak burst signals from 
Pop Ill/Pop II collapse. If, in the future, it will be possible 
to extract information on the Pop III GWB component, we 
can have a tool to investigate the formation history of these 
first stars. 
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